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Abstract 

An interpretation of scale-invariant multiplicity fluctuations inside hadronic 
jets is presented. It is based on the branching mechanism with the angular 
ordering of soft partons in sequential branchings. A relationship with fractal 
distributions is demonstrated. The model takes into account the finiteness of 
the number of particles produced in jets (finite energy) and leads to a good 
description of the multifractal fluctuations observed in e + e~ processes. 
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1 Introduction 



In the case of e + e~ annihilation processes the asymptotic collinear and infrared 
contributions to gluon cross sections can be described in Double Leading Log Ap- 
proximation (DLLA) by a Markov process (see Jl| for a review). This semi-classical 
description takes into account soft gluon interference effects on the basis of the 
angular ordering prescription when the parton emission is described by successive 
branchings and the available phase space is reduced to ever smaller angular regions 
(color coherence effects). The corresponding QCD master equation is an integral 
one and is based on Dokshitzer-Gribov-Lipatov-Altarelli-Parisi energy-distribution 
kernels. 

In the framework of this description, progress has been made in obtaining an- 
gular scale-invariantQ correlations between partons [0] (see || for a review). This 
approach, by conception, is a correlation one, based on the method of characteristic 
functionals. Hence, to derive directly measurable quantities such as normalized fac- 
torial moments or factorial cumulants, one needs to perform an integration of the 
correlation functions over the restricted phase-space region under study. This is pos- 
sible only after the use of many approximations and by identifying the phase-space 
regions which give the leading contributions ||. 

Apart from this problem, there are also more basic questions which restrict the 
direct comparison of the QCD correlation approach with experimental data. Firstly, 
the perturbative QCD calculations deal with an asymptotic behavior of the multipar- 
ton correlations valid only for very high energies. In an idealized jet, therefore, finite 
parton multiplicities in small phase-space bins and energy-momentum conservation 
effects are systematically ignored [0. This is one of the most important reasons 
leading to disagreements between the analytical predictions and e + e~ data @, |J. 
Secondly, the increase of the coupling constant for very small phase-space regions 
sets a limit for the validity of perturbative QCD. Thirdly, non-perturbative effects 
such as hadronization, resonance decays and Bose-Einstein correlations complicate 
the comparison of theoretical many-particle inclusive densities with the data even 
at LEP1 energies [f§§. 

In this paper, therefore, we propose a new way to study the correlations in terms 
of fluctuations in the multihadronic systems produced in high-energy processes. Be- 
ing based on a fluctuation approach to intermittency phenomenon (see recent re- 
views ||[7| on the subject of intermittency), the model a priori takes into account 
the finiteness of the number of particles in a single event (finite energies). In order 
to describe the local multiplicity fluctuations, we adopted the differential Markov 
equation for parton branching, which has been used to describe global multiplicity 
fluctuations in high-energy physics in |8HTT1 (see also references in 0). 

One of the key ideas of this approach is that, in contrast to a full phase space, 

lr The scale- invariance means that a dynamical characteristic X(l) of correlations/fluctuations at 
a given resolution I has the property X(Xl) = \~ L X(l) with a constant L characterizing dynamics 
of a multiparticle system. 
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a Markov branching process inside a small phase-space window of size 5 can be 
characterized by a probability P n (t, 5) of detecting n particles, in which a dependence 
on an evolution parameter t can be factorized from a phase-space ^-dependence 
(see Sect. |2|). As a consequence of this assumption, the scale- invariant fluctuations 
experimentally observed inside jets |L2| may be considered as a result of fractal 
phase-space distribution for each particle emitted in successive Markov branchings 
(Sect. |3|). Such an idea ultimately leads to the possibility of taking into account an 
inhomogeneity of the parton correlations inside a jet and a fairly good quantitative 



agreement with the e + e -annihilation data |12| and the JETSET 7.4 PS model |T3 
(Sect. D. 



2 Statistical treatment of branchings 
2.1 Global equation 

At high energies, gluons dominate the parton-parton cross section due to the large 
color factor and the infrared singularity. This means that a good high-energy ap- 
proximation should consider gluon branching only. For generality, however, we shall 
consider both gluons and quarks treating them as partons. 

Let t be the evolution parameter of the parton branching process. The t can 
be related to the parton virtuality Q and can be defined in the usual way |8|-i 1 



However, hereafter we shall never refer to the explicit form of this parameter and 
shall regard it as representing the extent of branching or just time. We assume that 
the branching process starts with t = and continues until some tmax determined 
by a QCD cut-off Qq. The initial condition for the probability distribution P n (t) of 
having n particles radiated by the initial one is 

P n=0 (t = 0) = 1, P n ^{t = 0) = 0. (1) 

In the following we shall see that, under the assumptions to be made below, the 
structure of local fluctuations depends neither on the particular definition of the 
evolution parameter, nor on the initial conditions. The purpose of the introduction 
of ([!]) is only to give an illustration of the notion of a typical initiation of the cascade 
and its further evolution. 

A probabilistic scheme |] of the perturbative parton shower is based on classical 
picture of the Markov chains of independent parton splittings. Each elementary 
parton decay depends on just the nearest "forefather". Let us define W\dt as the 
probability of branching a — > b + c during a small range of t, dt, according to one 
of the following decays: g — > gg, q — > qg, and g — > qq. The infinitesimal probability 
W\ in the leading log picture can be written as 



a.b 



Wi = E / az^p a ^ bc (z), (2) 
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where a s is the strong coupling constant and P a ^ c {z) are the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi energy-distribution kernels. The sum runs over all allowed 
parton branchings. For our simplified model, we will consider the case with a s = 
const, so that W\ is a (divergent) constant which does not depend on t. 

If there are n partons, the probability of the parton emission increases. Let W n d t 
be the probability that the parton system with multiplicity n radiates a new parton 
during the infinitesimal interval (t, t + At). Generally, W n depends on the parton 
multiplicity n. This can be taken into account as 

W n = w{n)W x , w(l) = 1, (3) 

where w(n) is a function of n reflecting an increase of the parton radiation. Then 
the Markov pure birth evolution equation for the multiplicity distribution P n (t) of 



having n partons at time t is well-known [|14 

BP (t) 

°-^p- = W n ^P n ^(t) - W n P n (t). (4) 

The solution of this equation is a global multiplicity distribution P n {t). Since the 
equation contains ingredients of perturbative QCD, an essential point is to regularize 
Wi and consider the branching evolution up to tmax determined by the QCD cut-off 
Qo. In order to compare the P n (^max) with the data, one usually resorts to the local 
parton-hadron duality hypothesis which states that n for partons is proportional to 
the n for observed hadrons. 

The differential equation (^) with constant (^-independent) vertex probabilities 
W n has been analyzed in [^[TT|]. One of the most popular solutions is a negative bi- 
nomial distribution which was derived in the leading log picture for gluons in quark 
jet |J. Deviations from this distribution observed in e + e~ annihilation data are 
usually connected with the shoulder structure and a quasi-oscillatory behavior of 
H q moments seen at Z° peak. Recently, however, the negative binomial distribution 
has been reestablished again: In was shown that the full-phase-space multiplicity 
distribution for e + e~ annihilation data can be well reproduced by a weighted su- 
perposition of two negative binomial distributions |15|, PI, associated to two- and 
multi-jet events or the contributions from bb and light flavored events. 

For the full phase space, there is no physical reason to define W n in momentum 
space: The global distribution is momentum independent. However, to obtain vari- 
ous momentum characteristics of particle spectra (such as the multiplicity of partons 
above a fixed momentum), a more complex integro-differential equations should be 
analyzed Jl], ||. Below we will discuss another way to include a momentum de- 
pendence using a statistical projection of equation (|]) into momentum phase-space 
domains. 

2.2 Local equation 

Obviously, if one counts only the particles produced within a certain small range of 
phase space, not all particles can be detected in it. Let 7 n (<5) be the probability of 
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observing one particle in a phase-space domain of size 5 if this particle belongs to 
the parton system of multiplicity n + 1 > 1 in the full phase space. We put 

< ln{8) < 1, (5) 

so that 

ln {5 = 0) = 0, 7 n(5 = fi) = 1, (6) 

where Q is the size of full phase space (8 < Q) which can be defined in 3-momentum 
phase space or, say, in rapidity, or azimuthal angle. 

For a phase-space element of size 8, if the system is in state n at time t, the 
probability of the transition n — > n + 1 in the interval (t, t + At) is 

ln (8)W n At + o(At,8), 

where, as before, W n describes the emission of one particle into the full phase- 
space Q and the factor j n (8) describes the probability of hitting 8 by this particle. 
The factorization property of the infinitesimal probability r y n {8)W n is an essential 
assumption used to simplify the structure of parton evolution. We also assume that 
the probability P n (t+ At, 8) of having n particles inside 8 at t+ At is fully determined 
by P n (t, 8) and P n ^i{t,8) in the same 8. In fact, for a particular ("angular") choice 
of phase space, this is consistent with the coherent branching with angular ordering, 
since the contribution of particles from phase-space regions outside of 8 is considered 
to be very small (see the discussion below). On the basis of these assumptions, one 
can write 

P n (t + At, 8) = 7n-i(«y)W n _iP n _!(t, 6) At + (1 - ln (5)W n At)P n (t, 6) + o(At, 5), 

where the second term is due to probability conservation. Then the corresponding 
Markov equation for the local multiplicity distribution P n (t, S) is 

BP (t S) 

] = Tn-iOW-iP^M) - ln {8) W n P n (t,5). (7) 

As we see, from the point of view of an observer counting particles in 5, the restriction 
of the phase-space domain looks as an effective suppression of the birth rate W n . 
(c.f.(^j)). Note also that, in contrast to ([|), equation (|7|) contains a momentum 
dependence via j n (8). 

It is necessary to note that condition (|5]) comes from a probabilistic interpretation 
of 7 n (5). Generally, as W n , this quantity can be larger than unity. However, if this 
is the case, we can carry out the following transition: 7 n (5) — > Zj n (8), where Z is a 
constant, so that the condition (§) for Zj n {8) can hold. As we shall see below, this 
regularization does not change the structure of observable fluctuations derived from 

(ED- 

Clearly, a possible non-linear nature of equation (0) renders its explicit solution 
very difficult. It can be solved in a straightforward manner only for some particular 
forms of the vertex probabilities W n and 7 n (5). 
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2.3 Phase-space property in the factorization scheme 

We will be interested in a general solution of (0) with respect to the possible behavior 
of the probability P n (t, 5) function of 7 n (#)- 
For n = 0, the solution can be easily obtained 

P (t, 5) = exp (- 70 (<5) J W d?j . (8) 

This exponential form of Pq is similar to the Sudakov form factor. In contrast to 
the full phase space, the integral contains the suppression factor 7 (5) taking into 
account the fact that a particle can be emitted outside of the small phase-space 
interval. 

The form of P n (t, S) for n > cannot be obtained without knowing the form of 
W n and j n (S). However, a phase-space structure of such a solution can be deduced 
in a general case. Since the basic idea of this approach is to factorize the phase- space 
and t-dependent component, let us look for the solution of (^) in the form 

P n (t, 5) = f n (t) Pn (5) + o(t, 5), n > 1, (9) 

where f n (t) is a 5- independent and p n (S) is t- independent well integrable functions. 
We assume that @ has a sense for any t at a sufficiently small 5. 
Using (H), (0) can be rewritten as 

M6) - *(*)&!, (10) 



Po(t,5) 



ww„_i(t) 



In-i(d) b n , n>2, (11) 
fo{t) = 1. (12) 



f n {t) + f n (t) W nln {5) ' 

Since we are looking for a solution at small 5, j n ($) has a small value. Therefore, 
b n can be approximated by the 5-independent constant, 

bn- f , {t) ■ (13) 

Further, the assumption (P) holds only if b n is independent of t for n > 2. For 
a given W n , (|13|) can be solved with respect to the form of f n {t)- However, the 
^-dependence of P n (S,t) has already been obtained. It reads 

P " iL()) -ln-x{5)b n 4^^1n-l{5)W n J-^, U > 1. (14) 



Pn-l(t,S) ~ ""/n-l(t) " '""^ y "~7;(*) : 

Let us remind that this relation is assumed to be possible only if S is small. In 
this case, the solution for P n (S, t) may be factorized as in (see an example in 



subsection f275[ ) 
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Of course, to study the distribution P n (t, S) as a function of 5 by means of 
factorial moments or cumulants might technically be a very difficult task. However, 
having in mind the bunching-parameter method , this distribution can easily 



be analyzed. Bunching parameters (BPs) r} q (t,S) are defined as 

They measure the deviation of the multiplicity distribution P n (t,5) from a Poisson 
one for which the BPs are equal to unity. Generally, in the case of no dynamical 
phase-space correlations, r) q (t,S) are 5-independent. 

The BP of an arbitrary order q for ([14]) can be written as 

<i,(t,s) = n,(t)v,(t), (is) 

w[q - 2) f q f q „i 
where rj q (5) depends only on the phase-space interval, 

9 7«-i(<*) / 10 x 



q- l7,_ 2 (5)" 



As we see, the structure of rj q (t, 6) is quite remarkable. It contains a t— dependent 
function r] q {t) constructed from unknown w(n) and f n (t), so that equation ([/]) itself 
can have strong non-linear property. However, since we study the fluctuations at 
ever smaller 5, this function is unrelevant: The property of the local fluctuations is 
fully determined by the ratio 7 g _i(5)/7 g __ 2 (5)- 

Note that while the original equation (^) is constructed from the divergent con- 
stants W n = w(n)Wi, the final result for the BPs does not contain them directly, 
since W\ cancels in (|17|). However, (|7|) contains them indirectly via f n (t). We can 
handle this problem since the regularization procedure 7 n (5) — > Z j n (5) discussed in 
subsection ^]2| does not change the BPs ( |TTD and, hence, the observable fluctuations. 
According to this, one can always redefine j n (6) as j n (S) — > H / 1 _1 7 n (5), so that W\ 
cancels already in (|^). 

2.4 Markov birth-death process 

The same phase-space behavior fllTf ) of the BPs can be obtained from a station- 
ary Markov birth-death evolution equation. For small 5, this process has to be 
characterized by the birth rate ^ n {5)W^ and the death rate W~ due to the fusion 
(absorption) processes such as gg — > g, qg — > q and qq — > g. These effects are not 
important for the full phase-space. However, for small 5, the values of 7 n (^)H / ^ t " and 
W~ can be comparable. The local equation reads 

8P (t 5) r i 

^ = 7n-l(5)Wn-lPn-l(t, 6) + W~ +1 P n+1 {t, 5) - [<y n (S) W+ + W~) P n (t, 6). 
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Assuming that for very small 8 the process is a stationary, dP n /dt ~ 0, one can 
derive (see details in lEOl) 



Pn(S) _ W+_l 



Pn-l(S) W~ 



ln-l(S), (19) 



which is similar to (0). Hence, BPs have the same form as (|16|), with the phase- 
space dependence as in flTgp. The only difference is that r) q (t) in flT7| ) does not 
depend on t and has the form: 



V*= T# • (20) 




Note that the stationary (equilibrium) regime is a strong assumption. It cannot 
be applied to the full phase space. For local phase-space domains, the physical 
situation is somewhat different: Each emitted parton increases the phase space for 
further emissions and the total phase space is expanded with increasing t. However, 
if one counts the particles inside a selected small phase-space window, one may 
assume that there is a little change in the density of partons inside S with increasing 
t and, hence, P n (t, 5) does not depend strongly on t. This assumption can be verified 
experimentally by observing t-independence of the BPs. 



2.5 Fully independent emission 

A simple example of the approach discussed above provides a fully independent 
particle emission. For this we should use the following assumptions: 

1) W n in (0) does not depend on n, i.e. w(n) = 1, W n = Wi, 

2) 7n(<5) does not depend on n, 7„(5) = 7(5). 

Under these conditions, equation ([7]) can be easily solved. The solution is a 
Poisson distribution, 

P n (5) = a n exp(-a)/n\, a = W 1 t 1 (5). (21) 

The behavior of this distribution at small 5 can be factorized as in @, 

P n (5)~(H/ 1 t) n 7 n W/n! + o( 7 n (5)), 

so that the corresponding BPs are unity. Note that for (pl|) this is true not only 
locally (5 — > 0), but also for any 5. For a uniform phase-space distribution, 7(5) is 
simply equal to 5/ A. 

Generally, an independent phase-space particle production can be characterized 
by any W n with 7„(<5) = 7(5). In this case the BPs are (^-independent constants. 
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3 Local fluctuations in the model 



3.1 Statistically averaged picture of a jet 

To study the phase-space dependence of fluctuations, the next step is to understand 
a possible behavior of y n (5) in 



We shall start our consideration with a simple two-dimensional model of a jet in 
angular intervals. Let us consider the first parton emitted at some angle with respect 
to the initial quark. Since we are interested in a picture averaged over all events, let 
Qq be the maximum possible size of solid angle, so that the first parton always has 
an angle inside the cone f2 (see Fig. p. After its emission, the first parton radiates 
the next one at some angle with respect to its own direction of flight. Generally, we 
assume that there is recoil effect and the first parton can change its direction after 
this radiation. In this case, the solid angular window available for both partons 
becomes larger and is equal to fix > Qq. The second parton then splits into two new 
partons at and so on. One can further simplify the model taking into the account 
angular ordering when available phase space is reduced for successive branchings. 
In this case fl ~ fix — fl2 — ■ ■ ■ fin- 
Let us tern to a more detailed description in one dimension. First, let us define 
as the polar angle between the directions of motion of the emitted and the parent 
parton. The single-particle distribution p(0) of the gluon bremsstrahlung can be 
approximated ||,j3| at small by 

p(<d)=C(Q ,a s )e-\ (22) 

integrating the overall distribution over the azimuthal angle around the quark direc- 
tion and momentum dependence. The ©-independent constant C(Q ,a s ) contains 
a transverse momentum cut-off Qo and a s which is treated here as a constant. The 
probability j (8Q) of finding the gluon inside the small interval (0 O — 50, ©o) near 
a jet opening angle ©o is 



7o(<50) oc / p(0)d© ~ 5Q V °, V = 1 (23) 
JBo-se 



for 50 — ^ 0. Note that this result does not depend essentially on the details of the 
density p(Q), since it has no singularity near O . We did not specify a coefficient of 
proportionality between 7 (<5O) and 5Q V °: As we have seen before, the phase-space 
dependence of the fluctuations does not depend on it. 

Now let us consider the behavior of 7i(0) for the second parton. Since we are 
interested in the probability of emission of this parton into (0 O — 5Q, O ) under the 
condition that the first parton is inside the same interval, there is a larger probability 
of hitting this interval by the second parton because of the collinear singularity. Now 
the major problems in the calculating 71 (0) are: 1) An ambiguity in the position 
of the first parton inside <50; 2) Singularity of p(0) near ~ gives a dominant 
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contribution. This leads to a very inhomogeneous phase-space distribution near 6 ; 
3) Requirement of the angular ordering. 

Due to the reasons quoted above, the calculation of 7 n (©) for n > 1 is even more 
difficult. We shall make no attempts to calculate 7„(G). In a general case, for small 
5, we assume 

j n (Se) oc 50 Vn , n>l, (24) 

where T> n are 50-independent constants controlling the collinear singularities to- 
gether with the angular ordering restrictions of the phase space available for particles 
on (n + l)th multiplicity stage. The latter effect decreases the available phase space 
for the next soft offspring partons that would increase the probability of detecting 
them inside 50. We assume, 

V >V l >V 2 >...>V oo . (25) 



In subsection [372] we shall give an interpretation of the behavior fl24 ) and (|25|) in 



terms of fractal distributions. Then we shall see that the behavior of 7 n (<50) for 
small 50 is the only simplest choice which allows to describe experimental data. In 
Sect. |] we shall proceed with the physical interpretation of these quantities. 

There are a number of special cases of interest: 

1) Monofractal fluctuations 

This case corresponds to the situation when the phase-space distributions for all 
cascade stages (except the initial one) have the same non-uniformity characterized 
by Pi, i.e., 

l0 (5O)oc5O v \ ln>0 (5O) oc 50 v \ (26) 
Making use of ([IB]) , the BPs are 

r] 2 (50) oc 50 Vl ~ Vo , V q >2(50) = const. (27) 



Hence, we obtain the monofractal behavior with 0I2 = T>o — T>\ |L7|, . 

For cascade branchings, such a situation can be considered as a highly unrealistic 
since it totally disregards that daughter partons have ever larger probability to be 
emitted inside 50 because of the correlations. Therefore, the monofractal type 
of intermittency possibly observed for some nucleus-nucleus reactions may mainly 



be attributed to other dynamical mechanisms plf , rather than to actual cascade 
processes with angular ordering. 

2) Multifractal fluctuations 

If particles on each cascade stage are distributed differently, then the cascade 
stage with the multiplicity n + 1 should be characterized by its own V n , i.e., 

7 „>o(^e) oc 50 v ". (28) 

The corresponding BPs are 

77,(50) oc 50~ a \ a q = V q _ 2 -V q _ x . (29) 
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An inverse relation for T> n reads 



V n = Vn 



n+l 

^2 on. 

i=2 



(30) 



According to [17, IS] , one has a multifractal behavior. An example of such a behavior 
will be given in subsection [3~4. 



3.2 Connection with fractals 

The simplicity of the model allows a natural connection of it with fractals. In this 
subsection we shall see that V n introduced in (|24D are nothing but fractal dimensions. 

First, let us remind a standard definition of a fractal distribution. Let us assume 
that there is a large number A^ ^ of particles distributed over a phase space with 
the topological (Euclidean) integer dimension D (D — 1,2,3). Let A/"(e) be the 
number of particles counted inside the phase-space domain with a linear size e. The 
number Af(e) and e are related as 

M{e) oc e°, e -> 0, (31) 

where T> is a fractal dimension, corresponding to the so-called box-counting (or mass, 
cluster etc.) dimension ||23|| . If the distribution is extremely inhomogeneous, T> has 
a non- integer value {T> < D). If particles were uniformly distributed over the phase 
space, T> is integer (T> = D). Therefore, T> is a very economical way to describe the 
extent of non-uniformity of a distribution near a given small phase-space region. 

It is easy to see that (|31|) also characterizes the probability p of observing one 
particle inside e: This probability is determined by the ratio of the number iV(e) of 
events of observing a particle inside e to the total number A^ t of events. Assuming 
that only one particle can be emitted in each event, one has 

Me) A/Ye) v . . 

P=TT i = T7 Li «^, e^O. (32) 

7V tot A/tot 

Now let us tern to the model. In fact, the 7 n (5B) has the same meaning as the 
p defined in (|32]). The index n in 7„(50) simply specifies the cascade stage n with 
the total n + l particles, so that D n stands the fractal dimension of the phase-space 
distribution of a single particle on each cascade stage. Then ( p3|) describes a uniform 
particle distribution near O (no collinear singularity!). For the second particle, 
there is no such a uniformity any more: The collinear singularity of the emission of 
the second particle is near 0o and this leads to a very inhomogeneous distribution 
in this region, so that 7i (59) oc 5Q Vl , where V x is a fractal dimension of this 
distribution (T>\ < T>q = 1) . For the next emissions, the distribution should be even 
more inhomogeneous since parent particles are already non-uniformly distributed 
due to the collinear singularities and the angular ordering. Finally this leads to the 



condition V n > V n+1 guessed in (|25|) 
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The T> n are the usual fractal dimensions. However, after many cascade steps with 
different T> n , one obtains a multifractal behavior fl29|) of the BPs. For a monofrac- 
tality (|2?D, the phase-space distribution for each particle in the cascade has to be 
characterized by a single fractal dimension for all n, T>q ^ T> = T>\ = T>2 = . . .. 



3.3 Connection with factorial- moment method 

A widely used means to study local fluctuations is based on the calculation of the 
normalized factorial moments 

el , An x (n(n-l)...(n-q + T)) 

F q (8G) = — , (33) 

where n is the number of particles inside a restricted phase-space interval 5Q, (...) 
is the average over all events. For non-statistical fluctuations, F q (5Q) depend on the 
size of the phase-space interval <50 as F q (5Q) ~ 5Q~^ q , where <ft q are intermittency 
indices. 

If the size of phase space is asymptotically small, then the following approximate 



relation between the FJ5Q) and the BPs holds [fT?, 18 



F q (5Q) ~ J] V q n +1 - n m- (34) 



n=2 



From (|34D and (2S), one has 

F q (6Q)cx6Q-^, <P Q = ^(g-n + l)^, (35) 

n=2 

or, taking into account the expression for at n , 

^= E(?-^+l)(^n-2-^n-l). (36) 
n=2 

The case of no dynamical correlation corresponds to <p q = 0. From fl36|) , it follows 
that the only possibility for this case is the condition 

V = 1 = V 1= V 2 = .... (37) 

i.e., the next emitted partons are distributed over available phase space purely ran- 
domly (uniformly). 

The model allows a simple way to connect the Renyi fractal dimension (see details 
in J7|) for factorial moments with the usual fractal dimensions T> q in our model. The 
Renyi fractal dimension D q is defined via (f> q , 

D q = D — -A^. (38) 
11 



From one has, 



£, = z>i-£ g n | 1 (^_ 2 - 2Vi), (39) 

where we take into account that the topological dimension .D is equal to T>q. From 
here one can again see that the monofractality (D q = const) is possible only if 
T> n _i = V n , for n > 1. A variation of D q with q for the multifractal case can be due 
to X> n _i ^ V n . 

In fact, the information about the fractal dimensions T> n can be extracted from 
the study of both D q (for factorial moments) or a q (for bunching parameters). 
However, the study of the BPs is the most direct way to obtain the information 
on V n : 

1) In contrast to the BPs, the power-like behavior of the normalized factorial 
moments holds only approximately for one dimensional variables because of a satu- 
ration effect for small rapidity intervals (see [fi,|7|,|l2"|). 

2) The BP 7] q of order q is a differential tool, resolving only the difference T> q _ 2 — 
T> q _i between the fractal dimensions D n (see (|29|)). In contrast, the normalized 
factorial moment F q of order q is an "integral" tool, which is sensitive to to all T> n 
with n — 1 < q. Because of the factor in the sum (p9|), the contribution from T> n at 
small n is the largest. Hence, small changes in the behavior of T> n for large n may 
be hidden due to contributions from T> n for small n. 



3.4 Experimental data 

The multifractal behavior (|29|) of BPs is characteristic for many different reactions 



ITf. For example, for rapidity variable with respect to the trust axis, BPs depend 



on the size of rapidity interval 5y as 

Vq (5y)=p' q 5y- a \ q > 2, (40) 

where and a q are positive constants. This can be considered as an evidence that 
local fluctuations have a scale-invariant structure, r] q (X5y) = X~ aq r] q (5y), i.e. the 
behavior is invariant under change of scale. 

Usually, the power law ( ^0|) is represented in terms of the number M = Y/Sy of 
bins of size 5y covering a full phase-space volume Y, so that ([40]) becomes 



n q {M)=(3' q M a ". (41) 

Taking the logarithm from both sides, the power law can be written as the linear 
expression 

ln Vq (M) = a q lnM + P q , P q = ln^. (42) 

For e + e~ annihilations, such a behavior has been observed for rapidity defined 
with respect to the thrust axis (see Fig. ||and |p!2| ,[17| , [l8| ). That the a q are not zero 



12 



and vary with q is a direct indication that the fluctuations in y are mult if r act al. 
Table [I] shows the values of ot q and j3 q obtained using a fit by fl4jp. To avoid trivial 
effects due to a bell-shaped structure of the multiplicity distribution at large M, the 
fit is limited to In M > 3 for q = 2 and In M > 2 for q > 2. 

Fig. H and Table [l] show the predictions of the JETSET 7.4 PS [D| model with 
the L3 default parameters The charged final-state hadrons were generated at 
91.2 GeV. The total number of events is 2.0M. The regions InM < 3 (for rj 2 and 
773) and InM < 2 (for 774) were excluded from the fits. Note also that x 2 test for 
the Monte Carlo is rather poor since, for the large statistics used, the behavior of 
r] q {M) shows a clear complex structure caused by the presence of resonance decay 
products and the points for different M are not statistically independent. 

Table [2] shows the fractal dimensions T> n obtained using The values of T* n 
decrease with increasing n, indicating that the degree of non-homogeneity of the 
distributions increases for particles emitted in the cascade later. 



4 Model predictions 

We have now set up a formalism that handles the local scale-invariant fluctuations 
inside a cascade. Qualitatively, the model proposed above can reproduce the power- 
like dependence of BPs observed in e + e~ data [T2j and other process [JT7] . 

A most direct prediction of this approach is that the power-like behavior of the 
BPs is energy independent: The local fluctuations are determined by 7 n (<5) in (|T5|). 
They, in turn, depend only on the fractal dimensions D n . As a result, parameters 
a q determining the phase-space fluctuations in ( |29"1) are t- independent. 

The model, however, has only low predictive power unless we reduce the number 
of free parameters V n in (|29"D . To do this, let us rewrite the V n as 

V n = V - A n , (43) 

so that positive A n represents the deviation of fluctuations from the trivial ones 
(A n = actually corresponds to the case of no correlations or uniform cascade 
distributions). We shall call the parameters A n as the strength of dynamical cor- 
relations on the n + 1 multiplicity stage of the branching. Since > T>„, we 
have 

A) = < A 1 < A 2 < . . . < A n . (44) 

The physical meaning of A n is rather clear: A n is determined by the collinear 
singularities of gluon emission and the extent of interference between soft partons 
leading to angular ordering. Generally, however, A n may absorb many other physi- 
cal effects in jet beyond DLLA. This quantity can incorporate effects from energy- 
momentum balance (recoil effect) in two-parton splittings, heavy quark production 
and non-perturbative effects: hadronization, resonance decay and Bose-Einstein cor- 
relations. Since contributions from these effects are poorly known and at present 
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cannot be taken into account in analytical calculations, below we shall make an 
attempt to treat A n on a general statistical ground. 

Several remarkable features of A n are immediately apparent: 

a) A n characterizes a single particle inside SO belonging to a system with n other 
particles already produced inside this interval at the previous cascade stages. 

b) Since A n is connected with correlations/fluctuations, one can consider it as a 
strength of "interaction" of a single particle with another. According to (f44"D, such 
an interaction becomes stronger with the increase of multiplicity n. 

These two features suggest that A n is analogous to the binding (pairing) energy 
per nucleon in nuclear physics. Using this analogy, the form of A n can be readily 
deduced without detailed information on correlations. 

Let us first consider the following two extreme cases: 

1) Since the Markov chain is based on two-particle splittings, one can assume 
that there exist positive correlations only between the particles a\ and a 2 of the 
two-particle splitting a\ — > a\ + a 2 , which is a basic element of the Markov chain. 
From a statistical point of view, the effect tends to make two partons more strongly 
bound in phase space, i.e., the probability that particles a± and a 2 occupy a very 
small phase-space bin is larger than that without dynamical correlations. After the 
next splitting of each particle, one has 2 two-particle pairs. For an (n + l)-particle 
system, the number of pairs stemming from the two-particle splittings is (n + l)/2, 
and we can write 

A n = A T ^-, (45) 

where A T is a constant describing the pair correlation in the case of two-particle 
correlations.^ Note that the applicability of fl4~5|) for odd n is only an approximation 
to make the correlations easy to handle. We shall correct this expression later. 
If only two-particle correlations ([45]) exist, then one obtains from (|43| ) and ( f2~9| ) 



T n + 1 



V n = V -A'^-, (46) 

a 2 = A T , a q > 3 = 0.5A T . (47) 
The behavior a q = 0.5a 2 has been found to correspond to multiplicity fluctuations 



in pp collisions \ T% . However, e + e~ data show a stronger multifractal signal. The 
behavior of (^) with v4 T =0.016 for the e + e~ data is shown in Fig. |] (A T > 0, A M = 
0). The value of A T is equal to a 2 taken from the experimental data (see Table [I]). 
The model fails to reproduce the n-dependence of T> n for data and JETSET model. 

2 The two-particle and multiparticle correlations introduced in our statistical model to describe 
the cascade have nothing to do with the two-particle and multiparticle correlations in the final- 
state hadrons measured by means of the two-particle and multiparticle correlation functions Q . 
We borrowed these terms following an analogy with the Weizsacker mass formula for the binding 
energy per nucleon in nuclear physics. 
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2) Now let us consider another limiting case of correlations. Let us assume that 
each particle of a given (n+l)th particle generation is attracted in equal extent by all 
of the other n particles already produced. There are exactly n(n + l)/2 interactions 
between n + 1 particles uniformly distributed over the small phase-space volume. 
(Such a uniformity must, of course, be treated as an average over all events.) Hence, 
the correlation strength is (see Fig. |5|) 

An = A u *P±H (48) 

where A M is a constant characterizing the correlation between any two particles. It 
completely determines many-particle correlations in such a system. 
Having made this simple assumption, one has 

V n = V -A M ^±^, (49) 

and, according to (|29|), the power-law indices for the BPs in the form 

a 2 = A M } a q > 3 = A M (q-l). (50) 

The result for A M = 0.016 is shown in Fig. | (A T = 0, A M > 0). As we see, this 
prediction is rather close to the experimental result. However, it still cannot give a 
satisfactory description of the data and JETSET model. In fact, such a disagreement 
is not a surprise since we systematically ignored the trivial fact that particles can 
interact with different strength. 

As was mentioned, to some extent, A n is analogous to the binding (pairing) 
energy per nucleon in nuclear physics. In fact, expression ( f4"5]) is analogous to 
the "volume" effect if the nuclear density is roughly constant. Then each nucleon 
has about the same number of neighbors and ( f4"5|) actually represents the short- 
range correlations. Then (f48|) is analogous to the Coulomb repulsion term in the 
Weizsacker mass formula which is proportional to |2H 



Z(Z-l) , N 

-a 2 (51) 

where Z is the number of protons and a = e 2 /4ir is the fine-structure constant 
of QED. The negative sign implies a reduction in binding energy. For QCD, of 
course, the Coulomb interaction is not the dominant part of the correlations and 
the introduced correlations should be attributed to other reasons. 

Following the same logic, A n can be constructed analogously to the semi-empirical 
Weizsacker mass formula by combining the different types of correlations and taking 
into account the obvious properties of the particle system in question. To see this, 
let us consider the following cascade chain: 

Oi —> («i + ai) —> (oi + 03) + «2 —> ■ ■ ■ , 
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where the a n represents a parton in independent sequential splittings. The parti- 
cles in parentheses are pairs arising due to two-particle splitting of parent particles 
on each stage. It is natural to assume that correlations between particles in the 
parentheses are different from those between the particles that have already been 
produced. For example, the particles in the pairs (ai,a 2 ) and (02,^3) produced on 
the three-particle stage can also be correlated, but to an extent different from those 
in the pair (a 1; a 3 ) which stem directly from the two-particle splitting. Thus to make 
a step towards a more realistic description, it is necessary to take into account a 
non-homogeneity of parton interactions in the cascade. 

First of all, let us describe the correlations between the particles in two-particle 
splittings. For this, we should take into account the odd-even effect in the two- 
particle correlations which is important for small n (this was dropped for simplicity 
in (f45"D). A corrected expression (|4^) reads as 

a t _ a t j (n+l)/2, for n = 1,3,5,... 

An ~ A X \ n/2, forn = 2,4,6,... [bZ) 

The next step is to take into account the multiparticle correlations arising be- 
tween the particles produced in the previous stages of the cascade. As before, to 
simplify our considerations, we assume that this kind of (multiparticle) correlations 
can be characterized by a single parameter A M responsible for the correlation be- 
tween any particles stemming from different parents. For any n-particle system, 
the form of these correlations can be obtained by subtracting from a term of the 
form (0), representing all possible pair correlations, a term like (|52"D describing 
two-particle correlations which are taken into account by (|52|). The final expression 
reads 

a m_ A M n ( n + 1 ) a m v / {n + l)/2, for n = 1,3,5,... 
A n= A 2 \ n/2, for n = 2,4,6,... {b6) 

The last step is to combine both contributions together, 

V n = V - Al - AX, (54) 

a 2 = A T , a 9 >3 = + Af_, - A^_ 2 - Af_ 2 . (55) 

Expressions ([52]), (|53"|), (|54|) and (|^) explicitly describe the behavior of the correla- 
tions in the cascade on the basis of the two parameters A T and A M . The parameter 
A T describes the correlation between particles stemming from the same parent parti- 
cle and A M characterizes the correlation between the particles coming from different 
parents. As in nuclear physics,^ we allow these constants to be adjustable and 
consider A T and A M as free parameters which can be evaluated from the fit. 

The parameters A M and A T can be obtained from the two experimental param- 
eters af X ^ and af X ^ describing the power-law behavior of BPs: 

A T = (56) 



3 In nuclear physics the situation is somewhat different: provides a "volume" binding effect 
with positive sign and has negative sign that implies a reduction in binding energy. 
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A M = «f P /2. 



(57) 



Further evolution of the V n and the a q can be predicted by the model according to 
fl54l) and (|55|) . For the e + e~ data presented in Table [I], one obtains A T = 0.016±0.004 
and A M = 0.021 ± 0.002. The predictions for V n are shown in Fig. | (A T , A M > 0). 
The dashed lines show the uncertainty in the behavior of T* n due to the statistical 
errors on A T and A M . Our predictions agree with the experimental data well. The 
agreement with the JETSET becomes better if one uses the values of and 0:3 
from the Monte-Carlo model to determine A T and A M . 

Note that expressions (p^) and (^) cannot be valid for asymptotically large n 
since the fractal dimensions T> n cannot be smaller then zero. 

5 Discussion of the model predictions 

One of the striking features of the results obtained is that good agreement between 
the model and the data is possible only if the value of A T is smaller than that of 
A M . This means that the binding effect between two particles from the same parent 
must be smaller than that between particles produced earlier from different parent 
particles, i.e., the particles originating from different parents have a larger chance 
of being emitted very close to each other. 

There are a number of possible explanations for this effect. If one believes that 
the model describes the perturbative QCD cascade, the reason for this may come 
directly from the color coherence effects. Indeed, the fact that A M > can be due 
to the angular ordering: For a given cascade stage with multiplicity n, collective 
correlation effects should exist between each particle due to the angular ordering 
history of the previous stages. Then the smallness of A T can be explained by recoil 
effects and the minimal value of the relative transverse momentum k±_ of decay 
products in the cascade evolution, in order to ensure that partons have enough time 
to radiate, in their turn, new offspring M. The latter effect leads to a restriction 
on the relative emission angels between the particles a\ and 02 in the two-particle 
splitting ai — > a\ + ci2- From a statistical point of view, the effect tends to make 
the two partons less tightly bound in phase space, i.e., the probability that both a\ 
and 02 particles occupy a very small phase-space bin is less than that without the 
restriction on the angle. If the reason for the condition A T < A M indeed comes from 
perturbative QCD, A^ has to be connected with the momentum transfer cut-off Qq 
that limits the relative k± and plays the role of an effective mass of a parton. 

On the other hand, it is reasonable to think that the proposed formulation of 
the branching process is sufficiently general and can utilize non-perturbative effects 
as well. In fact, the branching can be attributed to a certain degree to hadroniza- 
tion and resonance decay. Then, the multiparticle correlations can arise due to 
the color exchange between the partons at the end of the perturbative regime of 
QCD branching, necessary for parton discoloration. Furthermore, if the partons 
are replaced by hadrons, the large multiparticle correlations can be attributed to 
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Bose-Einstein interference between identical pions, since these particles are usually 
produced by different parent ones. Then the smallness of can be explained by 
an anti-correlation trend between decay products of resonances. 

Note also that the model can be used for various complex non-point-like pro- 
cesses. In this context, one can consider the evolution of the multiplicity distribu- 
tions for clusters, fireballs, resonances etc., taking into account peculiar features of 
these processes and introducing additional (or other) correlation terms in fl54|). 



6 Summary and conclusion 

In this paper we developed a new concept of local scale-invariant fluctuations in 
branching processes. In contrast to the approaches based on many-particle QCD cor- 
relation functions |2|,[j|] and phenomenological continuous densities []23j, we adopted 
a method based on single-particle probabilities (or single-particle probability densi- 
ties) for each cascade stage. They are characterized by the fractal dimensions T> n 
determining a non-uniformity in phase-space distributions for each particle emitted 
into a small phase-space domain. Such an idea simplifies the picture of phase-space 
organization of particles inside the cascade and allows us to take into account the 
finiteness of the number of particles in the cascade (or finite energy), QCD color 
coherence effects and a heterogeneity of correlations between partons belonging to 
the different cascade generations. 

The fractal dimensions T> n can be experimentally observed by calculating the 
BPs which resolve the difference £> n _i — T> n) according to (p9|). A less direct way to 
measure D n can be performed from the study of the normalized factorial moments 
(see©). 

The model suggests and makes experimentally accessible new physical quantities 
- pair correlation coefficients A M and A T determining T> n . The fact that none of 
these parameters are zero is due to the collinear singularities of the emission prob- 
abilities of soft partons. However, the way how these parameters determine the 
directly observable V n can be due to many reasons. In this paper we suggest such 
a relationship using a general statistical formalism, which, in terms of QCD, may 
absorb the details of coherence effects, high-order perturbative corrections, recoil 
effects and non-perturbative phenomena, i.e. all the effects which at present can be 
combined together only on the basis of Monte-Carlo simulations. We allow A M and 
A T to be adjustable that ultimately leads to good quantitative agreement with the 
local fluctuations in e + e~ processes. 

The model predicts that the experimentally observable parameters T> n determin- 
ing the scale- invariant behavior of BPs rj q (S) are energy independent. In addition, 
they do not depend on details of Markov equation in the full phase space. Both fea- 
tures follow from the factorization scheme used to derive the local fluctuations from 
the classical Markov branching equation for jet evolution and the angular ordering 
scheme which helps to construct the local version of this equation. Therefore, to 
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check this approach, precise data on the behavior of the BPs with energy are needed. 

Another model prediction is a suppression of positive correlations between the 
off-spring particles, A M > A T , a feature which can directly be detected from the 
study of g-dependence of the BPs. This prediction is also model dependent and the 
next step would be to understand how this effect can be changed if one uses another 
physical motivated parameterizations. 

In spite of its simplicity, the model describes the correlations between partons 
in branchings beyond the scope of the Leading Log Approximation of QCD. To 
leading order in InQ 2 , partons are free elementary quanta. Evidently, this situation 
corresponds to the particular case V n = 1 (for all n) in our scheme. Since the model 
is constructed on the basis of angular ordering, it takes advantage of the DLLA. 
However, for very small 5, the perturbative QCD ceases to be valid, since Q sets 
the limit of validity of the smallest bin size and perturbative expansion of QCD. 
Hence, dealing with very small phase-space intervals, our model goes beyond the 
perturbative QCD approximations studied in ||. At the same time, the model can 
take into account non-perturbative effects which are important if one goes beyond 
single-particle densities. It is evident that the price to pay for this progress in the 
description of multiparticle correlations inside jet is a purely statistical formalism 
eliminating the momentum dependence. 
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a q 


A, 


X7df 


a q 


P g 


X7df 




data 


JETSET 7.4 PS 


q 


= 2 


0.016 ±0.004 


0.244 ±0.018 


0.8/8 


0.0206 ±0.0005 


0.224 ±0.002 


2.4/11 


q 


= 3 


0.042 ± 0.003 


0.08 ±0.01 


8/12 


0.0434 ± 0.0007 


0.075 ±0.003 


22/13 


Q 


= 4 


0.062 ±0.004 


0.01 ±0.01 


9/12 


0.068 ±0.001 


-0.016 ±0.004 


36/10 


q 


= 5 


0.071 ±0.008 


-0.03 ±0.02 


14/11 


0.081 ±0.002 


-0.049 ± 0.004 


91/10 


q 


= 6 








0.072 ±0.002 


-0.019 ±0.005 


48/10 


q 


= 7 








0.088 ±0.003 


-0.053 ±0.006 


64/8 



Table 1: Fit results for r] q (M) obtained from the e + e data j7^/. The linear function 
/(Zjl) is used. 





data 


JETSET 7.4 PS 


n = 





1.0 


1.0 


n = 


1 


0.984 ±0.004 


0.9794 ±0.0005 


n = 


2 


0.942 ±0.005 


0.936 ±0.001 


n = 


3 


0.888 ±0.006 


0.868 ±0.001 


n = 


4 


0.81 ±0.01 


0.787 ±0.002 


n = 


5 




0.715 ±0.003 


n = 


6 




0.627 ±0.004 



Table 2: The values of fractal dimensions T> n obtained from the experimental data 
and JETSET 74 PS. (see §UJ and Table §}. 
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Figure 1: A schematic representation of the phase-space structure of branching inside 
jet. It makes use the angular ordering prescription: The structure of the cascade 
inside 5Q is determined by the "history" of this cascade inside the same 5Q. In- 
finitesimal probabilities W n (not shown) control the structure of the cascade for full 
phase space Q n . Local infinitesimal probabilities ^rWn determine the structure of 
cascade inside Sfl. 
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2345 2345 

In M In M 

Figure 2: BPs for rapidity defined with respect to the thrust axis for e + e~ processes. 
Here M = Y j by, "where Y is the size of full rapidity interval, Sy is the restricted 
rapidity interval. The data are reproduced from ftjjjjj . The lines represent the fit by 
ftm ) with the parameters presented in Table 
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1.1 




Figure 4: The behavior ofV n for e + e' - annihilation data, JETSET 7.4 PS and the 
model predictions for: a) Two-particle correlations (A T > 0,A M = 0); b) Mul- 
tiparticle correlations (A 1 = 0, A M > 0); c) Both two-particle and multiparticle 
correlations (A T ,A M > 0). 



n=l n=2 n=3 




Figure 5: A schematic representation of the multiparticle correlations for an (n + 1)- 
particle system (n = 1, 2, 3^. 
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